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Abstract 

cn ; 

I We investigate the trefoil and figure-eight knot universes from Bianchi-type I cosmology. 

f"*) . In particular, we construct several concrete models describing the knot universes related to 

£N| ' the cyclic universe and examine those cosmological features and properties in detail. Finally 

some examples of unknotted closed curves solutions (spiky and Mobius strip universes) are 

presented. 
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1 Introduction 



Inflation is one of the most important phenomena in modern cosmology and has been confirmed 
by recent observations on cosmic microwave background (CMB) radiation pQ. Furthermore, it 
is suggested by the cosmological and astronomical observations of Type la Supernovae [2], CMB 
radiation 1 , large scale structure (LSS) [3J, baryon acoustic oscillations (BAO) [3], and weak 
lensing [5] that the expansion of the current universe is accelerating. In order to explain the late 
time cosmic acceleration, we need to introduce so-called dark energy in the framework of general 
relativity or modify the gravitational theory, which can be regarded as a kind of geometrical dark 
energy (for reviews on dark energy, see, e.g., and for reviews on modified gravity, see, 



It is considered that there happened a Big Bang singularity in the early universe. In addition, 
at the dark energy dominated stage, the finite-time future singularities will occur [T5]-[23]. There 
also exists the possibility that a Big Crunch singularity will happen. To avoid such cosmological 
singularities, there are various proposals such as the cyclic universe [23]- [25] ( m other approach of 
the cyclic universe, see [27]), the ekpyrotic scenario [25], and the bouncing universe [25] . 

On the other hand, as a related theory to the cyclic universe, the trefoil and figure-eight knot 
universes have been explored in Ref. |30] . In the homogeneous and isotropic Friedmann-Lemaitre- 
Robertson- Walker (FLRW) and the homogeneous and anisotropic Bianchi-type I cosmologies, the 
geometrical description of these knot theories corresponds to oscillating solutions of the gravita- 
tional field equations. Note that the terms "the trefoil knot universe" and "the figure-eight knot 
universe" were introduced for the first time in Ref. [SO]- Moreover, the Weierstrass p(t), ((t) 
and a(t) functions and the Jacobian elliptic functions have been applied to solve several issues 
on astrophysics and cosmology 31 . In particular, very recently, by combining the reconstruction 
method in Refs. [T2J [22] [33J 133] with the Weierstrass and Jacobian elliptic functions, the equation 
of state (EoS) for the cyclic universes [35] and periodic generalizations of Chaplygin gas type mod- 
els [35]- [35] for dark energy [3j5] have been examined. This procedure can be considered to a novel 
approach to cosmological models in order to investigate the properties of dark energy. 

In this paper, we explore the cosmological features and properties of the trefoil and figure-eight 
knot universes from Bianchi-type I cosmology in detail. In particular, we construct several concrete 
models describing the trefoil and figure-eight knot universes based on Bianchi-type I spacetime. 
In our previous work |30) , the models of the knot universes from the homogeneous and isotropic 
FLRW spacetime were studied. By using the equivalent procedure, as continuous investigations, 
in this work we explicitly demonstrate that the knot universes can be constructed by Bianchi-type 
I spacetime. In other words, our purpose is to establish the formalism which can describe the knot 
universes. 

It is significant to emphasize that according to the recent cosmological data analysis pQ , it is 
implied that the universe is homogeneous and isotropic. In fact, however, recently the feature 
of anisotropy of cosmological phenomena such as anisotropic inflation [46] has also been studied 
in the literature. In such a cosmological sense, it can be regarded as reasonable to consider the 
anisotropic universe including Bianchi-type I spacetime. The units of the gravitational constant 
8ttG = c = 1 with G and c being the gravitational constant and the seed of light are used. 

The organization of the paper is as follows. In Sec. II, we explain the model and derive the basic 
equations. In Sec. Ill, we investigate the trefoil knot universe. Next, we study the figure-eight 
knot universe in Sec. IV. In Sec. V we present some unknotted closed curve solutions of the model. 
Finally, we give conclusions in Sec. VI. 



In this section we briefly review some basic facts about the Einstein's field equation. We start from 
the standard gravitational action (chosen units are c = 8itG = 1) 



e.g., na-nzi). 



2 The model 




(2.1) 
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where R is the Ricci scalar, A is the cosmological constant and L m is the matter Lagrangian. For 
a general metric g^, the line element is 

ds 2 = g^dx^dx", (/i,i/ = 0,l,2,3). (2.2) 

The corresponding Einstein field equations are given by 

R„u + (A - ^R)g^u = -k 2 T^, (2.3) 

where i? M „ is the Ricci tensor. This equation forms the mathematical basis of the theory of general 
relativity. In (I2.3[) . T^ v is the energy-momentum tensor of the matter field defined as 
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V = — (2-4) 



and satisfies the conservation equation 

V^T"" - 0, (2.5) 

where is the covariant derivative which is the relevant operator to smooth a tensor on a 
diffcrcntiable manifold. Ea. (l2.5[) yields the conservations of energy and momentums, corresponding 
to the independent variables involved. The general Einstein equation (|2.3p is a set of non-linear 
partial differential equations. We consider the Bianchi - I metric 

ds 2 = -dr 2 + A 2 dx\ + B 2 dx 2 2 + C 2 dxj, (2.6) 

where we assume that r = t/to, Xi = x[/xio, A, B, C are dimensionless (usually we put to = Xio = 
1). Here the metric potentials A, B and C are functions of r = t alone. This insures that the 
model is spatially homogeneous. The statistical volume for the anisotropic Bianchi type-I model 
can be written as 

V = ABC. (2.7) 

The Ricci scalar is 

n fA B C AB AC BC\ , nn . 

r = v jr - = 2 {a + b + c + ab + a^ + bc)> (2 ' 8) 

where A = dA/dr and so on. The non- vanishing components of Einstein tensor 

Gh = Rij - 0.5 gij R (2.9) 



are 



AB AC BC , oin , 
Go ° - AB + AC + BC> (2 - 10) 

Gaa = -^(i + S + fe) (2.11) 



Gbb = -Z 2 [a + C + 'ac)> (2 ' 12) 

«*> - - 2 (i+!+i)- 

We define a = (ABC)* as the average scale factor so that the average Hubble parameter may be 
defined as 
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^ = : = ; 7+d + t? ■ (2-14) 
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We write this average Hubble parameter H sometimes as 

H = ^(H 1 +H 2 + H 3 ), (2.15) 

where 

* = £ F 2 = §, F 3 = g (2.16) 

are the directional Hubble parameters in the directions of x\, x 2 and respectively. Hence we get 
the important relations 

A = A ef Hldt , B = B e$ H2d \ C = C ef Hadt , (2.17) 

where Aq 7 Bq,Cq are integration constants. The other important cosmological quantity is the 
deceleration parameter q, which for our model reads as 

(2-18) 

Next, we assume that the energy-momentum tensor of fluid has the form 

Tij = diag[To ,T 11 ,T 2 2,T 33 } = diag\p,-pi,-p 2 ,-pz]. (2.19) 

Here pi are the pressures along the Xi axes recpectively, p is the proper density of energy. Then 
the Einstein equations (with gravitational units, &ttG = 1 and c = 1) read as 

% " \R9ii = -Tij, (2.20) 
where we assumed A = 0. For the metric (|2.6p these equations take the form 

-P = 0, (2.21) 

Pi = 0, (2.22) 

P2 = 0, (2.23) 

P3 = 0. (2.24) 
In terms of the Hubble parameters this system takes the form 

H X H 2 + H 2 H 3 + H X H 3 - p = 0, (2.25) 

H 2 + H 3 +H$ + Hl+H 2 H 3 + pa = 0, (2.26) 

H i + H 1 +Hl + Hl + H 3 H 1 +p 2 = 0, (2.27) 

H 1 +H 2 +H 2 1 +Hj+H 1 H 2 +p 3 = 0. (2.28) 

Also we can introduce the three EoS parameters as 

Pi P2 P3 ,„ on , 

lji = — , lo 2 — — , lo 3 = — (2.29) 

P P P 

and the deceleration parameters 

AA BB CC , nnn . 

9l = -i7' q2 = ~^ Q3 = ~~fr- (2 - 30) 

Finally we want present the equation 

2H + 6H 2 = p-p, (2.31) 
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where 

Pi + P2 + pa , 

P = g (2-32) 

is the average pressure. Hence we can calculate the average papameter of the EoS as 

P + uj 2 + w 3 

w=- = . (2.33 

P 3 

Let us also we present the expression of R in terms of Hi . From (12.81) and (I2.16[) follows 

R = 2(H l +H 2 +H i + Hi + Hi + Hl + H X H 2 + H X H 3 + H 2 H 3 )) . (2.34) 
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Now we want present the knot and unknotted universe solutions of the system (|2.21l) - (|2.24l) or its 
equivalent (|2.25l) - (|2.28p . Consider some examples. 

3 The trefoil knot universe 

Our aim in this section is to construct simplest examples of the knot universes, namely, the trefoil 
knot universes. Consider examples. 

3.1 Example 1. 

Let us we assume that our universe is filled by the fluid with the following parametric EoS 

(3.1) 

(3.2) 

(3.3) 

(3.4) 
where 

Dx = (-12sin 2 (3r) + 36cos(3t) + 18cos 2 (3r))cos(2r) - 49sin(2r) 

(26/49 + cos(3r))sin(3r), (3.5) 

E 1 = sin(3r)(2 + cos(3r))sin(2r), (3.6) 
D 2 = -18sin(2r)cos 2 (3r) + (-49 sin(3r) cos(2r) - 36 sin(2r)) cos(3t) - 

26sin(3T)cos(2r) + 12 sin 2 (3r) sin(2r), (3.7) 

E 2 = sin(3r)(2 + cos(3r))cos(2r), (3.8) 
As = -30sin(3r)(2 + cos(3r))cos 2 (2r) -38sin(2r)(cos 2 (3T) - (27/38) sin 2 (3r) + 

(58/19) cos(3t) +40/19) cos(2r) + 30sin(3r) sin 2 (2r)(2 + cos(3t)), (3.9) 

E 3 = (2 + cos(3T)) 2 cos(2T)sin(2T), (3.10) 
Do = (6cos 2 (2t) - 6sin 2 (2r))cos 3 (3r) + (24cos 2 (2t) - 22 sin(3r) sin(2r) cos(2r) - 
24sin 2 (2T))cos 2 (3r) + ((-6sin 2 (3r) + 24) cos 2 (2t) - 52sin(3r) sin(2r) 
cos(2r) + (6sin 2 (3r) - 24) sin 2 (2r)) cos(3r) - (12(cos(2r) - (3/4) sin(3r) 

sin(2r)))(sin(3r) cos(2t) + (4/3) sin(2r)) sin(3r), (3.11) 

E = (2 + cos(3r)) 2 cos(2r)sin(2T)sin(3r). (3.12) 

Substituting these expressions for the pressures and density of energy into the system (|2.21l) - (|2.24p . 
we obtain the following its solution 

A = A + [2 + cos(3t)]cos(2t), (3.13) 
B = B + [2 + cos(3r)]sin(2r), (3.14) 
C = C* +sin(3T), (3.15) 
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where Aq,Bq,Co are some real constants. We see that this solution describes the trefoil knot. In 
fact the solution (|3.13[) - (I3.15I) is the parametric equation of the trefoil knot. In Fig. Q]we plot the 
trefoil knot for Eqs. (|3.13l) ~ (|3.15[) . where we assume 



A = B = Co = (3.16) 

and the initial conditions are A(Q) = 3, B(0) = C(0) = 0. The Hubble parameters for the solution 
|[3TT5 ]) -(13TT5 1) with (jXTCf read as 




Figure 1: The trefoil knot for Eqs. (pmj) - (PTT5j) . where A = B = C = 0. 

Hl = -2tan(2r)- 2sin(3T) (3.17) 
2 + cos(3t) 

H 2 = -2cot(2r)- 2sln(3T) (3.18) 
2 + cos(3t) 

H 3 = 3cot(3r). (3.19) 

In FigJ3]we plot the evolution of Hi for the solution p.l7[) ~ (|3.19l) with (I3.16p . It is interesting to 
study the evolution of the volume of the trefoil knot universe. For our case it is given by 

V = [2 + cos(3t)] 2 cos(2r) sin(2r) sin(3r). (3.20) 

In Fig. [3]we plot the evolution of the volume of the trefoil knot universe with respect to the cosmic 
time r for Eq. (|3.20[) with (|3.16[) . To get V > 0, we must consider Ao, Bo, Co > 0, if exactly e.g. 
as Aq > 3, Bo > 3, Co > 1- But below for simplicity we take the case (3.16). The other interesting 
quantity is the scalar curvature. For the trefoil knot solution (3.13)-(3.15), it has the form 

R = (6(12 sin 2 (3r) sin 2 (2r) - 28sin(3r) cos(2r) sin(2r) + 3sin 3 (3r) cos(2r) sin(2r) - 

- 12 sin 2 (3t) cos 2 (2t) - 8 cos(3r) sin 2 (2r) - 8 cos 2 (3t) sin 2 (2r) - 2 cos 3 (3r) sin 2 (2r) + 

+8 cos(3r) cos 2 (2r) + 8 cos 2 (3r) cos 2 (2t) + 2 cos 3 (3r) cos 2 (2r) - 

—52 sin(3r) cos(2r) cos(3r) sin(2r) — 19 cos(2r) sin(2r) sin(3r) cos 2 (3r) + 

+6 sin 2 (2r) sin 2 (3r) cos(3r) - 

-6 cos 2 (2r) sin 2 (3r) cos(3r)))/(sm(2r) cos(2r)(2 + cos(3r)) 2 sin(3r)). (3.21) 

In Fig|4] we plot the evolution of the R with respect of the cosmic time r. 

So we have shown that the universe can live in the trefoil knot orbit according to the solution 
(3.13)-(3.15). It is interesting to note that this trefoil knot solution admits infinite number acceler- 
ated and decelerated expansion phases of the universe. To show this, as an example let us consider 
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Figure 2: The evolution of the Hubble parameters for Eqs. (|3.17[) (|3.19|) . 

the solution for C from (3.13)-(3.15) that is C = Co + sin(3T). In this case we have C = — 9sui(3t) 
so that C > (accelerating phase) as r 6 (| + ^p, + ^ ZL ) and C < (decelerating phase) as 
r e (^f 2 -, § + ^f-) with the transion points C — 3cos(3t:;) = as 7$ = (0.57T + mr)/3, where rt is 
integer that is n = 0, ±1, ±2, ±3, .... 

3.2 Example 2. 

Now we consider the following parametric EoS 
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where 

Di = -sin 2 (2r)cos 2 (3r) + (-2cos(2t) - 4sm 2 (2-r) - 3 - sm(3r) sin(2r)) cos(3r) + 

+ sin(3r)sin(2r) - 4cos(2t) - sin 2 (3r) - 4sin 2 (2T), (3.26) 

Ei = 1, (3.27) 
D 2 = -(2 + cos 2 (3t))cos 2 (2t) - sin(3r)(-l + cos(3r)) cos(2r) + 

+ (-3 + 2sin(2r)) cos(3r) + 4sin(2r) - sin 2 (3r), (3.28) 

E 2 = 1, (3.29) 
D 3 = - (2 + cos 2 (3t)) cos 2 (2t) + (- sin(2r) cos 2 (3r) + (-4 sin(2r) - 2) cos(3t) + 3 sin(3T) - 

-4- 4sui(2t))cos(2t) + 2 sin(2-r) cos(3t) + (4 + 3 sin(3r)) sin(2r) - sin 2 (3r), (3.30) 

E 3 = 1, (3.31) 

D = (2 + cos(3T))(((2 + cos(3T))sin(2r)+sin(3T))cos(2r)+sin(3r)sin(2r)), (3.32) 

E Q = 1. (3.33) 

Substituting these expressions for the pressures and density of energy into the system (|3.22p - (|3.25l) . 
we obtain the following its solution 

Hi = [2 + cos(3r)] cos(2r) = 2 cos(2t) + 0.5[cos(5r) + cos(t)], (3.34) 
H 2 = [2 + cos(3t)] sin(2r) = 2sin(2r) + 0.5[sin(5r) - sin(r)], (3.35) 
H 3 = sin(3r). (3.36) 

We see that this solution again describes the trefoil knot but for the "coordinates" Hi. Note that 
the scale factors we can recovered from (|2.17l) . We get 

_ ^^gSin(2r)+0.1sin(5r)+0.5sin(T) (3 37) 

B = 5 oe -[ cos ( 2T )+ - lcos ( 5r )-°- 5cos ( r )] j (3.38) 

C = C Q e-? cos{3T \ (3.39) 

where Ao, Bq, Co are some real constants. In FigJS] we plot the evolution of A, B, C accordingly to 
(PO?)) - (PT31?)) and for the initial conditions A(0) = 1,-8(0) = e~ - 6 , C(0) = e" 1 / 3 , where we assume 
that Ao = Bq = Co = 1. For this example, the volume of the universe is given by 

V = y oe { sin ( 2r )+ - lsin ( 5r )+ - 5sin ( T )-[ cos ( 2T )+ - lcos ( 5T )-°' 5cos ( r )]-^ cos ( 3 ' r )} (3 40) 

The evolution of the volume for (|3.40[) is presented in Fig. [5] for Aq = Bq = Co = Vq = 1 and for 
the intial condition V(0) = e~ 14 / 15 . 
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400- 




Figure 4: The evolution of the R with respect of the cosmic time r for Eq. (|3,21[) 




Figure 5: The evolution of A, B, C accordingly to ([537)1 - ([339)1 . [t € 0..2tt]. 

The scalar curvature has the form 

R = (2 cos 2 (2t) + 2 sin 2 (2t) + 2 cos(2r) sin(2r)) cos 2 (3r) + (8 cos 2 (2r) + (2 sin(3r) + 4 + 
+8 sin(2r)) cos(2r) +6 + 8 sin 2 (2r) + (-4 + 2 sin(3r)) sin(2r)) cos(3r) + 8 cos 2 (2t) + 
+ (-2sin(3T) + 8 + 8sin(2T))cos(2r) + 8sin 2 (2r) + 

+ (-8 - 2sin(3<r)) sin(2r) + 2 sin 2 (3r). (3.41) 

In FigJTJ we plot the evolution of the R with respect of the cosmic time t. Finally we conclude 
that the Einstein equations for the Bianchi I type metric admit the trefoil knot solution of the form 
(3.34)-(3.36) or (3.37)-(3.39). These solutions describe the accelerated and decelerated phases of 
the expansion of the universe. 
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V 




Figure 6: The evolution of the volume for (|3.40p with Aq = Bq = Co = Vq = 1. 



3.3 Example 3. 

Now we present a new kind of the trefoil knot universes. Let the system (|2.2ip - (12.24p has the 
solution 

A = A + [2 + cn(3r)]cn(2r), (3.42) 

B = B + [2 + cn(3T)]sn(2r), (3.43) 

C = C + sn(3T), (3.44) 

where cn(<) = cn(i, k) and sn(i) = sn(i, k) are the Jacobian elliptic functions which are doubly 
periodic functions, k is the elliptic modulus. Fig. [8] shows the knotted closed curve corresponding 
to the solution (j3Ti2"|) - (f3~l4]) with (|3~TB)) . Substituting the formulas (|3T4^)) - (|3~i4]) into the system 
(|2.2ip - J2.24|) we get the corresponding expressions for p and pi that gives us the parametric EoS. 
This parametric EoS reads as 

Pi = ~, (3.45) 
P2 = ~, (3.46) 
P3 = (3.47) 
P = ^ (3-48) 
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Figure 7: The evolution of the R with respect of the cosmic time r for Eq. p.41j) 



where 
Dx = 



E-2 

D 3 



E?, 
Do 



En 



d d 
9fc 2 cn(3r, fc)sn 3 (3r, fc)sn(2r, fc) — 12— am(3r, fc)sn 2 (3r, fc)cn(2r, fc)— om(2r, k) — 

or or 

-4sn(2T, fc)((9/4)cn 3 (3r, k)k 2 + (9/2)cn 2 (3r, k)k 2 + 

+ ((45/4)^-am 2 (3r, fc) + cn 2 (2r, fc)fc 2 + |-am 2 (2T, fc))cn(3r, fc) + 2^-am 2 (2r, fc) + 

OT OT OT 







d 



(9/2)— am 2 (3T, fc) + 2cn 2 (2r, fc)fc 2 )sn(3r, fc) + 18cn(3r, fc)— am(3r, fc) 

OT 



dr 







(2 + cn(3r, fc))cn(2r, k)—am(2T, fc), 
£i = (2 + cn(3r, fc))sn(2r, fc)sn(3r,fc), 

d d 
Z?2 = 9fc 2 cn(3r, fc)sn 3 (3T, fc)cn(2T, fc) + 12— am(3T, fc)sn 2 (3r, k)—am(2T, fc)sn(2r, fc) - 

or or 

9 (9 

-9cn(2r, fc)cn 3 (3r, fc)fc 2 + 2cn 2 (3r, fc)fc 2 + ((4/9)— am 2 (2r, fc) + 5— am 2 (3r, fc) 



(3.49) 
(3.50) 
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t) 



-(4/9)fc 2 sn 2 (2r, fc))cn(3r, fc) + (8/9)— am 2 (2r, fc) + 2— am 2 (3r, fc) 



8t 



8t 



(3.51) 
(3.52) 



d 



-(8/9)fc 2 sn 2 (2r, fc))sn(3r, fc) - 

d d 
— 18cn(3r, k)— am(3r, fc)(2 + cn(3r, fc))— am(2r, fc)sn(2r, fc), 
or or 

sn(3r, fc)(2 + cn(3r, fc))cn(2r, fc), 

-4fc 2 sn(2r, fc)(2 + cti(3t, fc)) 2 cn 3 (2r, fc) - 30— am(2r, fc)sn(3r, fc)^-am(3r, fc) 

(2 + cn(3r, fc))cn 2 (2r, fc) + 4sn(2r, fc)(fc 2 (2 + cn 2 (3r, fc))sn 2 (2r, fc) + 

d d 
+ (-(9/2)— am 2 (3T, fc) + (9/2)fc 2 sn 2 (3r, fc) - 5— am 2 (2r, fc))cn 2 (3r, fc) + 

+ (-20— am 2 (2r, fc) + 9fc 2 sn 2 (3r, fc) - 9— am 2 (3r, fc))cn(3r, fc) + 

d d 
+ (27/4)sn 2 (3r, fc)— am(3r, fc) 2 - 20— am 2 (2r, fc))cn(2r, fc) + 

d d 
+30— am(3r, fc)sn(3T, fc)— am(2r, fc)sn 2 (2r, fc)(2 + cn(3r, fc)), 

(2 + cn(3r, fc)) 2 cn(2r, fc)sn(2r, fc), 

-4sn(3r, fc)sn(2r, fc)cn(2r, fc)(2 + cn(3r, fc)) 2 — am 2 (2r, fc) + 6— ara(r, fc)(2 + cn(3r, fc)) 

1 L 

(cn 2 (3r, fc) + 2cn(3r, fc) - sn 2 (3r, fc))(cn(2r, fc) - sn(2r, fc))(cn(2r, fc) + sn(2r, fc)) 
d_ 



(3.53) 
(3.54) 



am(2r, fc) - 18sn(3r, fc)sn(2r, fc)(-(l/2)sn 2 (3-r, fc) + cn 2 (3r, fc) 
d 



+2cn(3r, fc))cn(2r, fc)— am 2 (3r, fc), 

OT 

(2 + cn 2 (3r, fc))cn(2r, fc)sn(2r, fc)sn(3r, fc). 



(3.55) 
(3.56) 



Figure 8: The knotted closed curve corresponding to the solution (I3.42[) - (|3.44[) with p,16[) . t € 
[0,4tt], k = 1/3. 



The volume of the universe for the solution (|3.42l) - (|3.44j) with (|3 . 1 6[) looks like 

V = [2 + cn(3T)] 2 cn(2T)sn(2r)sn(3r). (3.57) 
The evolution of the volume for (|3.57[) is presented in Fig|H]The scalar curvature has the form 



3 -i 




-3 J 

Figure 9: The evolution of the volume of the trefoil knot universe with respect to the cosmic time 
r for Ea. (j3"37|) 
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R = (-8sn(2T,fc)sn(3r,fc)fc 2 (2 + cn(3r,fc)) 2 cn 3 (2T,fc) + 
+12dn(2r, fc)dn(3r, fc)(2 + cn(3r, fc))(cn 2 (3r, k) + 

+2cn(3T, k) - 3sn 2 (3r, fc))cn 2 (2r, k) - 18sn(3r, fc)(-(4/9)fc 2 (2 + cn(3r, fc)) 2 sn 2 (2r, k) + 

+ (-4fc 2 cn(3r, k) - 2cn 2 (3r, fc)fc 2 - dn 2 (3r, fc))sn 2 (3r, k) + 

+ (2 + cn(3r, fc))(cn 3 (3T, k)k 2 + 2cn 2 (3r, k)k 2 + 

+ (5dn 2 (3r, k) + (4/3)dn 2 (2r, fc))cn(3r, k) + 

+ (8/3)dn 2 (2r, k) + 2dn 2 (3r, fc)))sn(2r, fc)cn(2r, k) - 

-12dn(2r, fc)sn 2 (2r, fc)dn(3r, fc)(2 + cn(3r, fc))(cn 2 (3r, k) + 

+2cii(3t, k) - 3sn 2 (3r, fc)))/(cn(2r, fc)sn(2r, fc)(2 + cn(3r, fc)) 2 sn(3r, k)). (3.58) 
In FigJTU]we plot the evolution of the R with respect of the cosmic time r. 
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Figure 10: The evolution of the R with respect of the cosmic time r for Ea. (|3.58[) 



3.4 Example 4. 

Our fourth example is given by 

Hi = [2 + cn(3T)]cn(2T), (3.59) 

H 2 = [2 + cn(3T)]sn(2T), (3.60) 

H 3 = sn(3r) (3.61) 

which again the knotted closed curve in Fig. [3] but for the " coordinates" ifj . Note that the 
corresponding parametric EoS looks like 

Pi = (3.62) 

P2 = (3.63) 

P3 = ~, (3-64) 

^3 

P = % (3-65) 
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where 







Di = -(2 + cn(3r, fc)) 2 sn 2 (2r, fc) - sn(3r, fc)(-3— am(3r, k) + 2 + cn(3r, fc))sn(2r, fc) + 

or 

d d d 

+ (— 2cn(2r, fc)— am(2r, fc) — 3— am(3r, fc))cn(3r, fc) — 4cn(2r, fc) — am(2r, fc) — 
ar ar ar 



-sn 2 (3r, fc), 



£i = 1. 



(3.66) 
(3.67) 



d 



(2 + cn(3r, fc)) 2 cn 2 (2r, k) - sn(3r, fc)(-3— am(3r, fc) + 2 + cn(3r, fc))cn(2r, fc) + 



d d d 

+(— 3— am(3r, fc) + 2— am(2r, fc)sn(2r, fc))cn(3r, fc) + 4— am(2r, fc)sn(2r, fc) 

— sn 2 (3r, fc), 



£>3 



(3.68) 
(3.69) 



d 

-(2 + en(3r, fc)) 2 cn 2 (2r, fc) + (-sn(2r, fc)cn 2 (3r, fc) + (-4sra(2r, fc) - 2— am(2r, fc)) 

cn(3r, fc) + 3sn(3r, fc)— am(3r, fc) — 4— am(2r, fc) — 4sn(2r, fc))cn(2r, fc) + 

d d d 

+2— am(2r, fc)sn(2r, fc)cn(3r, fc) + (4— am(2r, fc) + 3sn(3r, fc)— am(3r, fc))sn(2r, fc) 
ar ar ar 



-sn 2 (3r, fc), 



1. 



(3.70) 
(3.71) 



£>o = (((2 + cn(3r, fc))sn(2r,fc) + sn(3r,fc))cn(2r,fc) + sn(2r,fc)sn(3T, fc)) 

(2 + cn(3r,fc)), 
E = 1. 



(3.72) 
(3.73) 



In Fig. QT]we plot the evolution of pi, p for (|3.62[) — (I3.65|i . The scalar curvature has the form 



Pi P2 Pi p 




Figure 11: The evolution of pi, p for (|3.62p - (|3.65l) . t £ [0, 2ir], k = 1/3, p(red), pi(blue), p 2 (green), 
Ps (black). 



R = 2(2 + cn(3r, fc)) 2 cn 2 (2r, fc) + (2sn(2r, fc)cn 2 (3r, fc) + (8sn(2r, fc) + 4dn(2r, fc) + 
+2sn(3r, fc))cn(3r, fc) + 8sn(2r, fc) + (4 - 6dn(3r, fc))sn(3r, fc) + 
+8dn(2r, fc))cn(2r, fc) + 2sn 2 (2r, fc)cn 2 (3r, fc) + (8sn 2 (2r, fc) + 
+ (-4dn(2r, fc) + 2sn(3r, fc))sn(2r, fc) + 6dn(3r, fc))cn(3r, fc) + 

+8sn 2 (2r, fc) + ((4 - 6dn(3r, fc))sn(3r, fc) - 8dn(2r, fc))sn(2r, fc) + 2sn(3r, fc) 2 (3.74) 
In FigfT2]we plot the evolution of the R with respect of the cosmic time r. 
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Figure 12: The evolution of the R with respect of the cosmic time r for Eq. (|3.74[) 

4 The figure-eight knot universe 

Our aim in this section is to demonstrate some examples the figure-eight knot universes for the 
Bianchi type I metric (|2.6I) . We give some particular figure-eight knot universe models. 

4.1 Example 1. 

Again, let us we assume that our universe is 

P 
Pi 

P2 
P3 



filled by the fluid with the following parametric EoS 

(4.1) 



E a ' 
_Ds 
E 9 ' 
Dip 
Eio ' 
Dii 
En ' 



(4.2) 
(4.3) 
(4.4) 
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where 



D 8 = (-2 sin(2r) cos(3r) - (3(2 + cos(2r))) sin(3r))(-2 sin(2r) sin(3r) + 
(3(2 + cos(2r)))cos(3T))sin(4r) + 12cos(3t)((2 + cos(2r)) cos(3t) - 
(2/3) sin(2r) sin(3r))(2 + cos(2r)) cos(4r) - (12(2 + cos(2r))) cos(4r) * 

((2 + cos(2r)) sin(3T) + (2/3) sin(2r) cos(3r)) sin(3r), (4.5) 

E s = (2 + cos(2T)) 2 cos(3T)sin(3T)sin(4r), (4.6) 
D 9 = ((72 + 36 cos(2t)) cos(4t) - 12 sin(4r) sin(2r)) cos(3t) - (29((24/29) * 

cos(4r) sin(2r) + (cos{2t) + 50/29) sin(4r))) sin(3r), (4.7) 

E 9 = sin(4r)(2 + cos(2r))sin(3r), (4.8) 
D 1Q = (-24cos(4T)sin(2r) + sin(4r)(-29 cos(2t) - 50))cos(3r) - 

(36((2 + cos(2r)) cos(4t) - (1/3) sin(4r) sin(2r))) sin(3r), (4.9) 

E 10 = sin(4r)(2 + cos(2T))cos(3r), (4.10) 
D n = -(30(2 + cos(2T)))sin(2r)cos(3r) 2 + sin(3r)(12 sin(2r) 2 - 196cos(2t) - 

180 - 53 cos(2r) 2 ) cos(3t) + 30sin(2r) sin(3r) 2 (2 + cos(2r)), (4.11) 

En = (2 + cos(2T)) 2 cos(3T)sin(3r). (4.12) 

Substituting these expressions for the pressuries and the density of energy into the system 
(|2~2l^ - ([2~2"4"|) . we obtain the following its solution [30] 



A = A + [2 + cos(2t)]cos(3t), 
B = B + [2 + cos(2r)]sin(3r), 
C = C* +sin(4T). 



(4.13) 
(4.14) 
(4.15) 



This solution is nothing but the parametric equation of the figure-eight knot as we can see from 
Fig. [T31 where we assume that Ao = Bq = Cq = and the initial conditions have the form 
A(0) = 3, B(0) = 0, C(0) = 0. And for that reason in [3D] we called such models as the figure-eight 
knot universes. Note that the "coordinates" A, B, C with (3.16) satisfy the equation 




Figure 13: The figure-eight knot for (j4TT3l) - (|IT31) with (|3~TFj) . 



4(/i - 2) 4 - 4(/i - 2) 2 + z 2 = 0, 



(4.16) 
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where h = 2 + cos(2r). Let us calculate the volume of the universe. For our case it is given by 

V = [2 + cos(2t)] 2 cos(3t) sin(3r) sin(4r), (4.17) 
where we used (|3.16[) . In Fig. 9 we present the evolution of the volume for the solution (|4.13[) - (|4.15l) 

V 

A A 




Figure 14: The evolution of the volume for the solution P~T5|| - P~T5|) with (l3"ToT) . t G [0, 2tt}. 

with (|3.16[) . The scalar curvature has the form 

R = ((24(cos(4r)cos(2r) - (3/2) sin(4r) sin(2r) + 2 cos(4r)))(2 + cos(2t)) cos 2 (3t) - 

-102sin(3T)(sin(4T) cos 2 (2t) + ((188/51) sin(4r) + (16/51) cos(4t) sin(2r)) cos(2t) + 

+ (32/51) cos(4r) sin(2r) + (172/51) sin(4r) - 

-(4/51) sin 2 (2r) sin(4r)) cos(3r) - 24sin 2 (3r)(cos(4r) cos(2r) - 

-(3/2)sm(4r) sin(2r) + 2cos(4t))(2 + cos(2r)))/(sin(3r) cos(3r) * 

*(2 + cos(2r)) 2 sin(4r)). (4.18) 
In FigfTS]we plot the evolution of the R with respect of the cosmic time r. So we found the figure- 



400 



200 



-200 




Figure 15: The evolution of the R with respect of the cosmic time r for Eq. (|4.18[) 

eight knot solution of the Einstein equations which again describe the accelerated and decelerated 
expansion phases of the universe. 
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4.2 Example 2 

Now we consider the system (j2.25p ~ (|2.28p . Let its solution is given by 

Hi = [2 + cos(2t)] cos(3t) = 2 cos(3t) + cos(5t) + cos(r), 
H 2 = [2 + cos(2r)] sin(3r) = 2 sin(3r) + sin(r) + sin(5r), 
H 3 = sin(4r). 

Then the coorresponding scale factors read as 

_ sin(3r)+0.2sin(5T)+sin(T) 

B = 5 oe -[f cos ( 3T )+ - 2cos ( 5T )+ cos (' r )] ; 

C* = C* oe -0.25cos(4r)^ 

For this solution the parametric EoS looks like 





D 


p = 


Eq' 




Di 


Pi = 


E x 




_D 2 


P2 = 


E 2 




D 3 


P3 = 


E 3 



where 



4.19) 
4.20) 
4.21) 



4.22) 
4.23) 
4.24) 



4.25) 
4.26) 
4.27) 
4.28) 



4.29) 
4.30) 

4.31) 
4.32) 

4.33) 
4.34) 



D a = (((2 + cos(2t)) sin(3r) + sin(4r)) cos(3t) + sin(3r) sin(4r))(2 + cos(2r)), 
E a = 1, 

Di = -(2 + cos(2r)) 2 sin 2 (3r) + (2 sin(2r) - 2sin(4r) - sin(4r) cos(2r)) sin(3r) - 

— 6cos(3t) — 3cos(3r) cos(2t) — 4cos(4t) — sin 2 (4r), 
Ei = 1, 

D 2 = -(2 + cos(2r)) 2 cos 2 (3r) + (2 sin(2r) - 2 sin(4r) - sin(4r) cos(2r)) cos(3r) - 

-4 cos(4t) + 6 sin(3r) + 3 sin(3r) cos(2r) - sin 2 (4t) , 
E 2 = 1, 

D 3 = -3sin(r) - 64 sin(r) cos 9 (r) + 36 sin(r) cos 5 (r) + 40 sin(r) cos 4 (r) + 4 sin(r) cos 3 (r) - 
-6sin(r) cos 2 (r) - 3sin(r) cos(r) - 25 cos 2 (r) + 5cos(r) - 40 cos 5 (r) - 64 cos 1 0(r) + 
+96 cos 8 (r) -84 cos 6 (r) +68 cos 4 (r) +26 cos 3 (r), (4.35) 

E 3 = 1. (4.36) 

In Fig. H6l we plot the EoS (|4.25|l - (|4.28p . For this example, the evolution of the volume of the 
universe is given by 

Y _ y^gf sin(3r)+0.2sin(5r)+sin(T)-| cos(3r)— 0.2cos(5r)-cos(r)— 0.25cos(4t) (4 37) 

The evolution of the volume is presented in FigfT7]for Aq = Bq = Co = Vq = 1 and for the intial 
condition V(Q) = e -127 / 60 . The scalar curvature has the form 

R = 2(2 + cos(2t)) 2 cos 2 (3t) + (2(2 + cos(2r)) 2 sin(3r) + (6 + 2 sin(4r)) cos(2r) + 
+ 12 - 4sin(2r) + 4sin(4r)) cos(3t) + 2(2 + cos(2t)) 2 sin 2 (3r) + 
+ ((-6 + 2 sin(4r)) cos(2t) + 4 sin(4r) - 4 sin(2r) - 12) * 

* sin(3r) + 2 sin 2 (4r) + 8 cos(4r) . (4.38) 

In Fig [18] we plot the evolution of the R with respect of the cosmic time r. Again we have 
shown that the Einstein equations admit the figure-eight knot solution and it again describe the 
accelerated and decelerated expansion phases of the universe. 



18 




4.3 Example 3. 

Now we present the figure-eight knot universe induced by the Jacobian elliptic functions. Let the 
system ([2T2"Tj) - ([2T2"I)) has the solution 

A = Ao + [2 + cn(2r)]cn(3r), (4.39) 
B = B + [2 + cn(2r)]sn(3r), (4.40) 
C = C + sn(4T). (4.41) 

Note that cn(£) and sn(t) are the doubly periodic Jacobian elliptic functions. FigfTH] shows the 
knotted closed curve corresponding to the solution (|4.39p ~ (l4.41l) with (|3.16[) . Substituting the 
formulas (|4.39p - (l4.41l) into the system (|2.2ip - (|2.24p we get the corresponding expressions for p 
and pi that gives us the parametric EoS. The evolution of the volume of the universe for (|3.16p 
reads as 

V = [2 + cn(2T)] 2 cn(3r)sn(3T)sn(4r). (4.42) 
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Figure 18: The evolution of the R with respect of the cosmic time r for Eq. (|4.38D 




Figure 19: The knotted closed curve corresponding to the solution (|4.39[) - (|4.41[) with (|3.16[l . 
t e [0,471-], k = 1/3. 



The scalar curvature has the form 

R = (-18sn(3r, fc)sn(4r, k)k 2 (2 + cn(2r, fc)) 2 cn 3 (3r, k) + (24(-(3/2)sn(2r, fc)dn(2r, fc)sn(4r, k) + 
+cn(4r, fc)dn(4r, fc)(2 + cn(2r, fc))))(2 + cn(2r, fc))dn(3r, fc)cn 2 (3r, k) - 
-(32(-(9/16)sn(4r, k)k 2 (2 + cn(2r, fc)) 2 sn 2 (3T, k) + ((cn 2 (4r, k)k 2 + (27/16)dn 2 (3r, k) + 
+dn 2 (4r, fc) - (l/2)/fc 2 sn 2 (2r, fc) + (l/2)dn 2 (2r, fc))cn 2 (2T, fc) + (-fc 2 sn 2 (2r, fc) + 
+ (27/4)dn 2 (3r, fc) + dn 2 (2r, fc) + 4dn 2 (4r, fc) + 4cn 2 (4r, fc)fc 2 )cn(2r, fc) + 4dn 2 (4r, fc) + 
+ (27/4)dn 2 (3r, fc) - (l/4)dn 2 (2r, fc)sn 2 (2r, fc) + 4cn 2 (4r, fc)fc 2 )sn(4r, fc) + cn(4r, fc)dn(4r, fc) 
dn(2r, fc)sn(2r, fc)(2 + cn(2r, fc))))sn(3r, fc)cn(3r, fc) - (24(-(3/2)sn(2r, fc)dn(2r, fc) 
sn(4r, fc) + cn(4r, fc)dn(4r, fc)(2 + cn(2r, fc))))(2 + cn(2r, fc))sn 2 (3r, fc)dn(3r, fc))/(cn(3r, fc) 
sn(3r, fc)(2 + cn(2r, fc)) 2 sn(4r, fc)). (4.43) 

In Figl2(I]we plot the evolution of the R with respect of the cosmic time r. 
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Figure 20: The evolution of the R with respect of the cosmic time r for Eq. (l4.43|) 



4.4 Example 4. 

We now consider the following solution of the system (|2.25l) - (|2.28[) : 

Hi = [2 + cn(2T)]cn(3r), 
H 2 = [2 + cn(2r)]sn(3r), 
H 3 = sn(4r) 



(4.44) 
(4.45) 
(4.46) 



which again the trefoil knot universe as shown in FigQ2J] but for the "coordinates" Hi. The 
corresponding parametric EoS reads as 



P 
Pi 

P2 
P3 



Do 

Eo' 
Di 
Ei' 

_Lh 
E2 ' 

E3 ' 



(4.47) 
(4.48) 
(4.49) 
(4.50) 
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where 
Do 

E Q 



D 2 



(((2 + cn(2r, fc))sn(3r, fc) + sn(4r, k))cn(3r, k) + sn(3r, k)sn(4r, fc)) 

(2 + cn(2T, k)), 

1, 

9 9 
2— am(2r, fc)sn(2r, fc)sn(3r, fc) — (3(2 + en(2r, fc)))cn(3r, fc)— am(3r, fc) 

-4cn(4r, fc)— am(4r, fc) - (2 + cn(2r, fc)) 2 sn(3T, fc) 2 - sn(4r, fc) 2 - 
ar 

-(2 + ctj(2t, fc))sn(3r, fc)sn(4r, fc), 



(4.51) 
(4.52) 



1, 



(4.53) 
(4.54) 



9 9 

-4ctiUt, k)—am(4T, k) + 2— am(2r, k)sn(2r, fc)cn(3r, fc) + (3(2 + cn(2r, fc))) 

OT OT 



0_ 



am(3r, fc)sn(3r, fc) - sn 2 (4r, fc) - (2 + cn(2r, fc)) 2 cn 2 (3r, fc) 



-(2 + cn(2r, fc))cn(3r, fc)sn(4r, fc), 
£ 2 = 1, 

D 3 = -(2 + cn(2r, fc)) 2 cn 2 (3r, fc) + (-sn(3r, fc)cn 2 (2r, fc) + (-4sn(3r,fc) - 

9 9 9 

— 3— am(3r, fc))cn(2r, fc) + 2— am(2r, fc)sn(2r, fc) — 6— am(3r, fc) 
or 9t ot 

9 

— 4sn(3r, fc))cn(3r, fc) + 3— am(3r, fc)sn(3r, k)cn(2r, fc) + 

9 9 
+ (6— am(3r, fc) + 2— am(2r, k)sn(2r, fc))sn(3r, fc) — sn 2 (4r, fc), 

OT OT 

E 3 = 1. 
Its plot we give in Fig|2"Tl 



(4.55) 
(4.56) 



(4.57) 
(4.58) 




Figure 21: The plot of the EoS (|Q71) - (|43H| . i € [0,2tt], fc = 1/3, p(red), pi(blue), p 2 (green), 
P3 (black). 

The scalar curvature has the form 

R = 2(2 + cn(2r, fc)) 2 cn(3r, fc) 2 + (2(2 + cn(2r, fc)) 2 sn(3r, fc) + (6dn(3r, fc) + 2sn(4r, fc))cn(2r, fc) + 
+12dn(3r, fc) - 4dn(2r, fc)sn(2r, fc) + 4sn(4r, fc))cn(3r, fc) + 2(2 + cn(2r, fc)) 2 sn(3r, fc) 2 + 
+((-6dn(3r, fc) + 2sn(4r, fc))cn(2r, fc) + 4sn(4r, fc) - 4dn(2r, fc)sn(2r, fc) - 12dn(3r, fc)) * 
*sn(3r, fc) + 2sn(4r, fc) 2 + 8cn(4r, fc)dn(4r, fc). (4.59) 

In Figl22]we plot the evolution of the R with respect of the cosmic time r. 



22 



R 



-20- 



40- 



20- 



30- 



10- 



10 - 








Figure 22: The evolution of the R with respect of the cosmic time r for Eq. (|4.59D 

5 Other unknotted models of the universe 

In this section we would like to present some unknotted but closed curve solutions of the Einstein 
equation for the Bianchi I type metric. As an examples we consider the spiky and Mobious strip 
universe solutions. 

5.1 Spiky universe solutions 

Our aim in this subsection is to present some unknotted closed curve solutions namely the spiky 
universe solutions. 

5.1.1 Example 1 

Let our universe is filled by the fluid with the following parametric EoS 



P2 = 



P3 = 



Pi 



P 



Dip 
Eio 
Dn 
En 



(5.4) 



(5.2) 



(5.3) 



(5.1) 
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where 

As = 



^8 

D 9 

E 9 
Dio 



— [asin((n — l)r)(n — 1) + a(n — 1) sin(r)] [a cos((n — l)r)(n — 1) — 

a(n — 1) cos(r)] sin(r) + [a cos((n — l)r)(n — 1) — a(n — 1) cos(r)] cos(r) x 

a cos((n — l)r) + a(n — 1) cos(r)] — [a sin((n — 1)t)(71 — 1) + a(n — 1) sin(r)] x 

asin((n — l)r) — a(n — 1) sin(r)] cos(r), 

a cos((n — l)r) + a(n — 1) cos(r)][a sin((n — l)r) — a(n — 1) sin(r)] sin(r), 



^10 



J5n 



(5.5) 
(5.6) 



—a sin((n — l)r)(n — l) 2 + a(n — 1) sin(r)] sin(r) — [a sin((n — l)r) — a(n — 1) sin(r)] sin(r) + 
acos((n — l)r)(n — 1) — a(n — 1) cos(t)] cos(r), (5.7) 
a sin((n — l)r) — a(n — 1) sin(r)] sin(r), (5-8) 
a cos((n — 1)t) + a(n — 1) cos(t)] sin(r) + sin(r)[a cos((n — l)r)(n — l) 2 + a(n — 1) cos(r)] + 



\Y - a(n - 1) cos(r)] + 



a sin((n — l)r)(n — 1) + a(n — 1) sin(r)] cos(r), 
-[acos((n — 1)t) + a(n — 1) cos(r)] sin(r), 
asin((n — l)r) — a(n — 1) sin(r)][— a cos((n — l)r)(n — l) 2 
a cos((n — l)r) + a(n — 1) cos(r)][— a sin((n — l)r)(n — l) 2 + a(n — 1) sin(r)] — 
a sin((n — l)r)(n — 1) + a(n — 1) sin(r)] [a cos((n — l)r)(n — 1) — a(rt — 1) cos(t)], 
a cos((rt — 1)t) + a(n — 1) cos(r)][a sin((rt — 1)t) — a(n — 1) sin(r)]. 

Substituting these expressions for the pressuries and the density of energy into the system (|2.21[1 - 
(|2.24l) . we obtain the following its solution 

A = acos[(n — l)r] + a(n — 1) cos[r], (5.13) 
B = asin[(n - l)r] - a(n — l)sin[r], (5-14) 
C = sin(r). (5.15) 

It is the spiky like solution so that such solutions we call the spike universe. Its plot presented in 
Fig(53]for the initial conditions A(0) — an = 10, B(0) = 0, C(0) = 0. Let us calculate the volume 




(5.16) 



(5.9) 
(5.10) 



(5.11) 
(5.12) 



Figure 23: The spiky universe for (I5.13[) - (I5.15[) . n = 10, a = 1. 

of this universe. It is given by 

V = a 2 [cos[(n — 1)t] + (n — 1) cos[r]] [sin[(n — l)r] — (n — 1) sin[r]] sin(r). 

In FielMl shown the evolution of the volume for (|5.16p . n = 10, a = 1. The scalar curvature has 
the form 

R = (-2cos(r)(n - l)cos 2 ((n- 1)t) + ((6(4/3 - 2n + n 2 )) sin(r) sin((n - 1)t) - 

-(2((n - 2) cos(t) 2 + sin 2 ( 7")(n 2 — in + 4)))(n- l))cos((n- l)r) + 2cos(r)(sin 2 ((n - 1)t) 4 
+ sin(r)(n 2 - 4n + 6) sin((n - l)r) + (cos(r) 2 - 5sin 2 (r))(n - l))(n - l))/((cos((n - l)r) + 



+ cos(r)(n — 1))(— sin((n — l)r) + (n — 1) sin(r)) sin(r)). 



(5.17) 
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Figure 24: The evolution of the volume for (5.16), n = 10, a = 1. 
In Figl^5] we plot the evolution of the R with respect of the cosmic time r. In this example, we 
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Figure 25: The evolution of the R with respect of the cosmic time r for Eq. (|5.17|) 

have shown that the Einstein equations admit the spike-like solution. We can show this solution 
describes the accelerated and decelerated expansion phases of the universe. 

5.1.2 Example 2 

The system (|2.25l) - (|2.28|) admits the following solution 

Hi = a cos[(n — l)r] + a(n — 1) cos[r], (5.18) 
i?2 = asin[(n — l)r] — a(n — 1) sin[r], (5.19) 
H 3 = sin(r). (5.20) 
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The corresponding EoS takes the form 





■Pia 


p = 






Dl3 


Pi = 


E13 




Du 


P2 = 


Eu 




Du 


P3 = 


E15 



(5.21) 
(5.22) 
(5.23) 
(5.24) 



where 



D12 = [a cos((n — 1)t) + a(n — 1) cos(r)][asin((n — l)r) + [1 — a(n — 1)] sin(r)] + 

[a sin((n — l)r) — a(n — 1) sin(r)] sin(r), (5.25) 

E12 = 1, (5.26) 

-D13 = — l)[cos((n — l)r) — cos(r)] + cos(t) + [asin((n — 1)t) — a(n — 1) sin(r)] 2 + 

[a sin((n - l)r) + [1 - a(n - 1)] sin(r)] sin(r), (5.27) 

E 13 = 1, (5.28) 

Z?i4 = — asin((n — l)r)(n — 1) — a(n — 1) sin(r) + cos(t) + [acos((n — l)r) + a(n — 1) cos(r)] 2 + 

sin(T) 2 + [acos((n — l)r) + a * (n — 1) cos(t)] sin(r), (5.29) 

#14 = 1, (5.30) 

D15 — a(n — l)[cos((n — l)r) — cos(t) — sin((n — l)r) — sin(r)] + 

[a sin((n — l)r) — a(n — 1) sin(r)] 2 + [a cos((n — 1)t) + a(n — 1) cos(r)] 2 + 

[a cos((n — 1)t) + a(n — 1) cos(r)][a sin((n — 1)t) — a(n — 1) sin(r)], (5.31) 

E 15 = 1. (5.32) 

The scalar curvature has the form 

R = 2a 2 cos((n - l)r) 2 + 4a((l/2)asin((n - l)r) + (1/2 + (-(l/2)n + l/2)a) sin(r) + 

+(n- l)(acos(r) + l/2))cos((n - l)r) + 2a 2 sin((n - l)r) 2 + 2a((l + (2 - 2n)a) sin(r) + 
+ (acos(T) - l)(n - l))sin((n - l)r) + (2 + 2a 2 (n - l) 2 + (2 - 2n)a) sin(r) 2 - (2(n - l))a * 
*(1 + (-1 + a(n - 1)) cos(r)) sin(r) + (2(a 2 (n - l) 2 cos(r) + 1 + a(-n + 1))) cos(r). (5.33) 

In Fig 1211 we plot the evolution of the R with respect of the cosmic time r. 
5.1.3 Example 3 

Our next solution for the system (|2.25[) - (|2.28[) is given by 

Hi = a cos[(n — l)r] — a(n — 1) cos[r], (5.34) 
H2 = asin[(n — l)r] — a(n — 1) sin[r], (5.35) 
H 3 = sin(r). (5.36) 

In Fig. 27 we plot this spiky type solution. The corresponding EoS takes the form 

(5.37) 
(5.38) 
(5.39) 
(5.40) 
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Eie 




Dn 
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En 
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Eig 
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Figure 26: The evolution of the R with respect of the cosmic time r for Eq. (|5.33D 




Figure 27: The evolution of the spiky type solution (5.34)-(5.36) with n = 10, a = 1. 

where 

Z?i6 = [a cos((n — l)r) — a(n — 1) cos(r)][asin((n — 1)t) + [1 — a(n — 1)] sin(r)] + 

[a sin((n — l)r) — a(n — 1) sin(r)] sin(r), (5-41) 

Eta = 1, (5.42) 

Z?i7 = a(n — l)[cos((n — l)r) — cos(t)] + cos(t) + [a sin((n — l)r) — a(n — 1) sin(r)] 2 + 

[a sin((n - l)r) + [1 - a(n - 1)] sin(-r)] sin(r), (5.43) 

Etr = 1, (5.44) 

£>i8 = — asin((n — l)r)(n — 1) + a(n — 1) sin(r) + cos(t) + [acos((n — l)r) — a(n — 1) cos(r)] 2 + 
sin(r) 2 + [a cos((n — l)r) — a * (n — 1) cos(r)] sin(r), (5.45) 

E 1S = 1, (5.46) 

Z?i9 = a(n — l)[cos((n — l)r) — cos(t) — sin((n — l)r) + sin(r)] + 

[a sin((n — l)r) — a(n — 1) sin(r)] 2 + [a cos((rt — 1)t) — a(rt — 1) cos(r)] 2 + 

[acos((n — 1)t) — a(n — 1) cos(r)][asm((n — l)r) — a(n — 1) sin(r)], (5-47) 

£ 19 = 1. (5.48) 
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The scalar curvature has the form 



R = 2a 2 cos((n- 1)t) 2 - (4(-(l/2)asin((ra - l)r) + (-1/2 + ((l/2)n - l/2)a) sin(r) + 
+(n- l)(-l/2 + acos(r))))o;cos((n - l)r) + 2a 2 sin((n - l)r) 2 - 
-(2((-l + (-2 + 2n)a) sin(r) + (a cos(r) + l)(n - l)))a sin((ra - 1)t) + 
+(2 + 2a 2 (n - l) 2 + (-2n + 2)a) sin(r) 2 + (2(n - 1))(1 + (-1 + a(n - 1)) * 
* cos(r))a sin(r) + 2 cos(r )(a 2 (n - l) 2 cos(r) + 1 + (1 - n)a). (5.49) 

In FigJ2H]we plot the evolution of the R with respect of the cosmic time r. 

300 -i 
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Figure 28: The evolution of the R with respect of the cosmic time r for Eq. (|5.49|) 



5.2 Mobius strip universe solutions 

If we consider the model with the " cosmological constant", then the systems (|2.2ip - (|2.24p and 
([2T2l)]) - ([2T38")) take the form, respectively 

-p-A = 0, (5.50) 

•Pi -A = 0, (5.51) 

p 2 -A = 0, (5.52) 

p 3 -A = (5.53) 
and 

H X H 2 + H 2 H 3 + H X H Z - p - A = 0, (5.54) 

H 2 + H 3 + Hi + Hl + H 2 H3+ Pl -A = 0, (5.55) 

H 3 + H 1 +Hl + Hf + H 3 H 1 +p 2 - A = 0, (5.56) 

Hx+H 2 + Hl + Hl + HxH 2 +p 3 - A = 0. (5.57) 

Now we want to present some solutions of these systems. Consider examples. 



AB 
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CA 


AB A 


BC 


f CA 
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CA 
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' CA ' 
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AB 
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5.2.1 Example 1 

One of the simplest solutions of (|5.50p ~ (|5.53p is given by 



.4 
B 



An + ( 1 H A COS — I COS T 



Bn - { 1 + ^Acos^ ) sill r. 



1 T 

C = Co + -Asm-. 



(5.58) 
(5.59) 
(5.60) 



It is the parametric equation of the Mobius strip and, hence, such model we call the Mobius strip 
universe. Its plot was presented in Fig |29l The evolution of the volume of the Mobius strip universe 




Figure 29: The plot of the Mobius strip universe for (|5.58[) - (15.60[) with (|3.16l) and t = — > 2ir and 
A =[-1.1] 



for (|535|l - (|5T50l with (|3~Td| reads as 



1 t\ t 
V = 0.5A I 1 H — A cos— cos r sin r sin — . 

2 2/ 2 



The evolution of the volume with (l3.16[) and a = A = 1 presented in FigMfl 
The corresponding EoS takes the form 



P 
Pi 

1>2 
P3 



D20 

E20 ' 
D21 
E21 ' 

E22 ' 
D23 
E23 ' 



(5.61) 



(5.62) 
(5.63) 
(5.64) 
(5.65) 
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Figure 30: The evolution of the volume of the Mobius strip universe for (|5.58l) - (|5.60p with (13. 16)) 
and a = A = 1. 



where 

D 20 = [~A sin(~) cos(r) + (1 + cos(~)) sin(r)] [ j A sin(~) sin(r) - (1 + ~A cos(^) cos(r)] x 
[C + ~Asin(~)] + ^[-iAsin(I)sin(T) + (1 + ^Acos(^)) cos(r)] cos(^) x 
[An + (1 + i)Acos(^)cos(r)] + ^[-^Asin(^) cos(r) - (1 + ^Acos(^)) sin(r)] cos(^) x 
[B + (1 + - )A cos(^)) sin(r)] - A[A + (1 + i)A cos(^)) cos(r)] x 

[Bo + (1 + iAcos(^))sin(r)][C + ^sin(^)], (5.66) 
E 20 = [A + (l + iAcos(^))cos(r)][S + (l + iAco S (^))sin(r)][C + iAsin(^)], (5.67) 
D21 = [Co + iAsin(^)][-isin(r)Acos(^) - ^A sin(^) cos(r) - (1 + ^Acos(^)) sin(r)] - 
- [Bo + (1 + ~ A cos(^)) sin(r)] sin(^) + - [- jA sin(^) sin(r) + (1 + i A cos(|)) x 
cos(r)] cos(^) - k[Bo + (1 + ^Acos(^)) S in(r)][C + ^Asin(^)], (5.68) 
E21 = [Bo + (l + iAcos(^))sin(r)][C + iAsin(l)], (5.69) 

D22 = ~[Ao + (1 + ~ cos(~)) cos(r)]Asin(^) + [C + ~Asin(~)][-^ cos(r) cos(^) + 

1 t It 1 1 t 

-Asin(-)sin(r) - (1 + -Acos(-)) cos(r)] + -[- jAsin(-) cos(r) - 

(1 + lAcos(^))sin(r)]Acos(^) - A[A Q + (1 + ^A cos(^)) cos(r)] [C + \a sin(J)], (5.70) 

£22 = [A + (l + iAcos(^))cos(r)][Co + iAsin(^)], (5.71) 

^23 = [Bo + (1 + ^Acos(^))sin(r)][-- cos(r) cos(^) + ^Asin(^) sin(r) - (1 + ^Acos(^)) x 

It A t 1 r 

cos(t)] + [A + (1 + -Acos(-))cos(r)][--sin(r) cos(-) - -Asin(-) cos(r) - 

(1 + ^Acos(^)) sin(r)] - [-Asin(^) cos(r) + (1 + ^Acos(^)) sin(r)] x 

[-^Asin(I) sin(r) + (1 + ^Acos(^)) cos(r)] - 

A[A a + (1 + iAcos(^)) cos(r)][B + (1 + ^Acos(^)) sin(r)], (5.72) 
2 2 30 2 2 

E23 = [B + (l + iAcos(^))sin(r)][A + (l + iAcos(^))cos(r)]. (5.73) 



The scalar curvature has the form 

R = ((-2sin(r) 2 A 2 + 2cos(t) 2 A 2 )cos((1/2)t) 3 + 

+(-8 sin(r) 2 A - 17 cos(r) sin(r) sin((l/2)r) A 2 + 

+8cos(r) 2 A) cos((1/2)t) 2 + ((6 sin(r) 2 A 2 - 6 cos(r) 2 A 2 ) sin((l/2)T) 2 - 

-60cos(r)sin(r)sin((l/2)r)A + 8cos(r) 2 - 8 sin(-r) 2 ) cos((l/2)r) + 

+ (sin((l/2)r) 2 A 2 cos(r) sin(r) + (12 sin(r) 2 A - 12 cos(r) 2 A) sin((l/2)r) - 

-52 cos(r) sin(r)) sm((l/2)r))/(sm(r) cos(t)(2 + Acos((1/2)t)) 2 sin((l/2)r)). (5.74) 

In FigEUwe plot the evolution of the R with respect of the cosmic time r. In this subsubsection, 




Figure 31: The evolution of the R with respect of the cosmic time r for Ea. (|5.74|) 

we have shown that the Einstein equations have the Mobius strip universe solution. Again we can 
show that this solution describes the accelerated and decelerated expansion phases of the universe. 

5.2.2 Example 2 

For the system (|5.54[) — f|5.57|) the Mobious solution reads as 



Hi = (l 


+ \ Acos (^ cos(t), 


(5.75) 


Hi = (l 


+ \ AC0S (^ sin ( T )> 


(5.76) 


H 3 = -A 

J 2 


sin( — ). 


(5.77) 


The corresponding EoS takes the form 






P 


£24 

E24. 


(5.78) 


Pi 


D25 
E25 ' 


(5.79) 


P2 


J>26 
E26 ' 


(5.80) 


P3 


D27 

E27 ' 


(5.81) 
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where 



£>24 

E24 

D 26 
D 27 



[1 + iAcos(^)] 2 cos(t) sin(r) + -[1 + - cos(^)][sin(r) + cos(r)] sin(^) - A, (5.82 



= 1, 



-Asin(^)sin(r) + [1 + iAcos(^)][cos(r) + -sin(r) sin(^)] + ^Acos(^) + 



[l + iAcos(^)] 2 sin 2 (r) + ^sin 2 (^)-A, 



= 1, 



(5.83 

- 

(5.84 
(5.85 



-^Asin(^)cos(r) - [1 + ^A cos(^)] [sin(r) + | cos(r) sin(J)] + ^Acos(^) + 



A. 



[l + iAcos(^)] 2 cos 2 (r) + ^sin 2 (^) 
1, 

A T A T 

I 1 + 2 cos ( 2 > ~ 1 sin ( 2 ^ ^ cos ^ + sin ( T )l 

A T 

[1+2 cos ( 2 )] 2 1 1 + C0S ( T ) sin ( T )l - A > 



E27 — 1. 
The scalar curvature has the form 

Jj> _ /1 /n\4 2/-„„„2/_\ i „ ; ,,-J 



(l/2)A 2 (cos 2 (r) + sin^r) + cos(r) sin(r)) cos 2 ((l/2)r) + 
+(l/2((cos(r) + sin (r))Asin((l/2)r) + 4cos 2 (r) + (2 + 4sin(r)) cos(r) + 
+1 - 2sin(r) +4sin 2 (r)))Acos((l/2)T) + (1/2)A 2 sin 2 ((l/2)r) + 
+(l/2(cos(r) + sm(r)))Asm((l/2)r) + 

+2cos 2 (r) + (1/2(4 + 4 sin(r))) cos(r) - 2 sin(r) + 2sin 2 (r) 



(5.86 
(5.87 

(5.88 
(5.89 



(5.90) 



In Fig|32]we plot the evolution of the R with respect of the cosmic time r. 




Figure 32: The evolution of the R with respect of the cosmic time r for Eq. l|5.90|) 
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5.3 Other examples of Mobius strip like universes induced by Jacobian 
elliptic functions 

5.3.1 Example 1 

Now we want present some solutions in terms of the Jacobian elliptic functions. In fact, the system 
(|5.50p - (|5.53[) has the following particular solution 



B 



C = 



A 




(5.91) 



(5.92) 



(5.93) 



The corresponding EoS takes the form 



Pi 



P 



D 2S 
E28 ' 
D29 
E29 
D30 
E30 
D31 
E31 



(5.94) 



(5.95) 



(5.96) 



(5.97) 
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where 

^28 = 



E28 
D 29 



E29 
D30 



E30 



E: 
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(5.98) 
(5.99) 
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-Adn — sn Acn — sn— 1 [-Adn— sn— snT — (1 H — Acn—) cnT dnr] x 



4 L 4 

cn^ dn^ - A[B Q + (1 + ^Acn^) siit][C + 5 Ami], 
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= -L4 + (1 + ^A cnl) cnr^A dn 2 I sn- + ^A cn 2 I snj 
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It It 
-A cn-) cnT][C + -Asn-], 
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2 2 8 2 2 
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—A dn — cn— cnT 
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-A dn— sn— dnr snT + (1 + -A cnI)(sn 2 T - dnV) cnr] + [A + (1 + - A cn-) cnT] x 
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-A cn— snrisn — 
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dnV) 
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1 A T \ 2 1 

—A cn— cn r snT 
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1 T T 

A dn— sn— cnT dnT — (1 
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1 t t It 
-A dn— sn— cnT + (1 H — A cn— ) dnT snrl x 
4 2 2 1 2 2' J 

1 t t It It 

[-A dn— sn— snT — (H — A cn— ) cnT dnT] — A[t4 + (1 H — A cn— ) cnT] x 
4 2 2 2 2 2 2 

[Bo + (1 + ^A cnl) snT], 

It It 

[B + (1 + 2 A cn 2) snr Po + (1 + 2 A cn 2) cnr ]' 



(5.104) 
(5.105) 



The evolution of the volume of the universe for (|3.16p reads as {Aq = Bo = Co = 0) 

2 



1 



1 A T 

-Acn— 

2 2 



cnT snT sn— . 

2 



(5.106) 



The evolution of the volume with (13. 16)) and A = 1 presented in Figji 
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Figure 33: The evolution of the volume of the trefoil knot universe with respect to the cosmic time 



r for Ea. (j5TU5j) 

The scalar curvature has the form 

R = (2A(2fc 2 cn((l/2)r,fc) + Afc 2 cn 2 ((l/2)r,fc) + 

+ (l/2)Adn 2 ((l/2)T, fc))cn(r, fc)sn(r, fc)sn((l/2)r, fc) 3 - 
-6Adn(r, fc)dn((l/2)r, fc)(cn(r, fc) - sn(r, fc))(cn(r, fc) + 
+sn(T, fc))(2 + Acn((l/2)r, fc))sn 2 ((l/2)r, fc) - 
-(4(2 + Acn((l/2)r, fc)))((l/4)cn 3 ((l/2)T, fc)fc 2 A + 

+ (l/2)cn 2 ((l/2)T, fc)fc 2 + A(3dn 2 (r, fc) + cn 2 (r, fc)fc 2 + (5/4)dn 2 ((l/2)r, fc) - 
-sn 2 (r, fc)fc 2 )cn((l/2)r, fc) + 6dn 2 (r, fc) + (l/2)dn 2 ((l/2)r, fc) - 
-2sn 2 (r, fc)fc 2 + 2cn 2 (r, fc)fc 2 )cn(r, fc)sn(r, fc)sn((l/2)r, fc) + 
+2cn((l/2)r, fc)dn((l/2)r, fc)dn(r, fc)(cn(r, fc) - sn(r, fc))(cn(r, fc) + 
+sn(r, fc))(2 + Acn((l/2)r, fc)) 2 )/(cn(r, fc)sn(r, fc) * 



In Fig|M]we plot the evolution of the R with respect of the cosmic time r. 
5.3.2 Example 2 

Similarly, we can show that the system (|5.54j) - (|5.57[) has the following solution 



*(2 + Acn 2 



((l/2)r,fc))sn((l/2)r,fc)). 



(5.107) 




(5.108) 



(5.109) 



(5.110) 
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Figure 34: The evolution of the R with respect of the cosmic time r for Eq. (|5.107[) 



The corresponding EoS takes the form 

P 
Pi 

P2 
P3 

where 



D32 

E32 
D33 
E33 
D34 
E34 ' 
D35 
E35 ' 



It A 1 t t 

D32 = [1 + - A cn-] 2 cnr snr + — [1 + -A cn-][snr + cnr] sn- - A, 

E32 — 1, 

1 T T 1 T ./Y T I T 

D 33 = -Adn- sn-[l - snr] + [1 + -A cn-][cnr dnr + — snr sn-] + [1 + -Acn-] 2 sn 

IaW^-A, 
4 2 

^33 = 1, 

1 t t It A tItt 

D34 = A dn— sn— cnr — (1 H — A cn— )[dnr snT H cnr sn— 1 H — A cn— dn— 

4 2 2 K 2 2 n 2 2 J 4 2 2 

[1 + cn J 2 cn 2 T + - A 2 sn 2 ^ - A, 
E 3 i = 1 , 

1 t t It 
B35 = — —A dn— sn— [cnT + snr] + [1 + — A cn— ][cnT — snTjdnT + 

1 T 

[1 + - A cn-] 2 [sn 2 T + cn 2 T + cnT snT] - A, 



E35 — 1- 



5.111) 
5.112) 
5.113) 
5.114) 



5.115) 
5.116) 



V 



5.117) 
5.118) 

5.119) 
5.120) 

5.121) 
5.122) 
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The scalar curvature has the form 

R = (l/2)A 2 (cn 2 (r,fc) + sn 2 (r,fc) + cn(r, fc)sn(r, fc))cn 2 ((l/2)r, fc) + 

+ (l/2(A(sn(r, fc) + cn(r, fc))sn((l/2)r, fc) + dn((l/2)r, k) + 4cn 2 (r, k) + 

+ (4sn(r, k) + 2dn(r, fc))cn(r, k) + 4sn 2 (r, k) - 

-2dn(r, fc)sn(r, fc)))Acn((l/2)r, k) + (l/2)A 2 sn 2 ((l/2)r, k) - 

-(l/2)A(dn((l/2)r, fc) - 2)(sn(r, k) + 

+cn(r, fc))sn((l/2)r, fc) + 2cn 2 (r, fc) + (l/2(4dn(r, fc) + 

+4sn(r, fc)))cn(r, fc) - 2sn(r, fc)(-sn(r, fc) + dn(r, fc)). (5.123) 
In Fig J3S] we plot the evolution of the R with respect of the cosmic time r. 
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Figure 35: The evolution of the R with respect of the cosmic time r for Eg. (|5. 1241) 



5.4 Integrable models 

The system (2.21)-(2.24) contents 4 equations for 7 unknown functions. This means we can add 3 
new additional equations. It gives us in particular to construct integrable Bianchi models. In this 
subsection we present two examples of such integrable models. 

5.4.1 Euler top equation 

Let us we assume that A = Ai, B = A-2,C = A3 obey the Euler top equation. The simple Euler 
top equation reads as 

ii = A 2 A 3 , (5.124) 

A 2 = A 3 A U (5.125) 

A 3 = A X A 2 . (5.126) 

The system of equations (2.21)-(2.24) and (5.124)-(5.126) we call the Bianchi I-Euler model which 
is integrable. 

5.4.2 Heisenberg ferromagnet equation 

Our second example is the Heisenberg ferromagnet equation (HFE). Here we assume that the 
variables A = S\,B = S 2 ,C = S3 satisfy the equation iS t + ^[S, W] = 0, iW x + u[S, W] = 0, 



37 



where S = SjCj, W — WiUi and Ui are Pauli's matrices. It the principal chiral type equation with 
the U — —i\S, V = — un^yj W.. Note this principal chiral type equation is the particular case 
of the following (2+l)-dimensional M-XCIX equation [37] (see also [IB] ) 



(that equivalent) 



iSt + 0.5([S,S y ]+uS) x + ~[S,W] = 0, 

UJ 

u x - 0.55 ■ [S x , S y ] = 0, 

iW x + u[S,W] = 



iS t + 0.5[S,S xy ] + uS x + -[S,W] = 0, 



UJ 

u x - 0.55 ■ [5c, Sy] = 0, 
iW x +w[S,W] = 0. 

This M-XCIX equation is integrable by the following Lax representation [37] 

$ x = 17$, 

$i = 2A$ y + 

with (5 2 = /) 

C/ = -iA5, 

V = AVi + —^—W - -W 
A + uj uj 



and 



5 



W 



s 3 s- 
s+ -s 3/ 

Vi = 0.25([5,5 B ] +w5), 
\W+ -W 3 



5.127) 

5.128) 
5.129) 



5.130) 

5.131) 
5.132) 



5.133) 
5.134) 



5.135) 
5.136) 

5.137) 
5.138) 
5.139) 



6 Conclusion 



In the present paper, we have constructed several concrete models describing the trefoil and figure- 
eight knot universes from Bianchi-type I cosmology and examined the cosmological features and 
properties in detail. 

To realize the cyclic universes, it is necessary to a non-canonical scalar field with non well-defined 
vacuum in the context of the quantum field theory or extend gravity, e.g., with adding higher order 
derivative terms and f(R) gravity [25]. Indeed, however, these modified gravity theories have to 
satisfy the tests on the solar system scale as well as cosmological constraints so that those can be 
alternative gravitational theories to general relativity. The significant cosmological consequence 
of our approach is that we have shown the possibility to obtain the knot universes related to the 
cyclic universes from Bianchi-type I spacetime within general relativity. 

Furthermore, recently it has been pointed out that the asymmetry of the EoS for the universe 
can lead to cosmological hysteresis [26]. On the other hand, Bianchi-type I spacetime describes 
the spatially anisotropic cosmology and hence the EoS for the universe has the asymmetry in the 
oscillating process through the expanding and contracting behaviors. Accordingly, it is considered 
that in the constructed models of the knot universes cosmological hysteresis could occur. The 
observation of this phenomenon in our models is one of our future works on the knot universes. 

Finally, it should be remarked that by summarizing the results of our previous [30] and this 
works, the knot universes describing the cyclic universes can be realized from the homogeneous and 
isotropic FLRW spacetime as well as the homogeneous and anisotropic Bianchi-type I cosmology. 
In these series of works, the formulations of model construction method of the knot universes have 
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been established. Thus, it can be expected that the presented formalism is useful to realize the 
universes with other features from both the isotropic and anisotropic spacetimes. 

Finally we would like to note that all solutions presented above describe the accelerated and 
decelerated expansion phases of the universe. 
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